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ABSTRACT 

We present a polarimetric map of a 20' x 20' area toward the Galactic center. The polarization 
of point sources has been measured in the J, and Ks bands using the near-infrared polarimet- 
ric camera SIRPOL on the 1.4 m telescope IRSF. One percent or better accuracy of polarization 
degree is achieved for sources with J < 14.5, H < 13.5, and Ks < 12.0. Comparing the Stokes 
parameters between high extinction stars and relatively low extinction ones, we have obtained a 
polarization originating from magnetically aligned dust grains at the central region of our Galaxy 
of at most 1—2 kpc. The distribution of the position angles shows a peak at ~20°, nearly parallel 
to the Galactic plane, suggesting a toroidal magnetic configuration. The derived direction of the 
magnetic field is in good agreement with that obtained from far-infrared/submillimeter observa- 
tions, which detect polarized thermal emission from dust in the molecular clouds at the Galactic 
center. Our results show that by subtracting foreground components, near-infrared polarimetry 
allows investigation of the magnetic field structure at the Galactic center. 

Subject headings: Galaxyicenter — dust, extinction — polarization — infrared:ISM — ISM:magnetic fields 



1. INTRODUCTION 

The discovery of polarized radio emission ex- 
tending perpendicular to the G alactic plane, such 



as no nthermal radio filament s (lYusef-Zadeh et aL. 
1984 ) and polarized plumes ( Tsuboi et aLlllQSGj T 
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has provided early evidence for a substantial 
poloidal component of the magnetic field in 
the central region of our Galaxy. High resolu- 
tion observations at radio wavelengths have re- 
vealed more filaments of similar orientation s (e.g., 
LaRosa et al.ll2QQ4tlYusef-Zadeh et"aIll2QQ4l ) . where 



the magnetic fields were found to be predomi- 
nantly aligned along the filaments. The accumu- 
lation of these observational results has led to the 
hypothesis that most of the volume of the Galactic 
center (GC) is permeated by a poloidal magnetic 
field. 

Evidence for the existence of a toroidal mag- 
netic structure comes from far-infrared (FIR) and 
submillimeter (sub- mm) observations, which de- 
tect polarized thermal emission from magnetically 
aligned dust grains. The rotation axis of the dust 
grain aligns with the magnetic field, and thus the 
polarization of the thermal emission indicates the 
magnetic field direction (the measured direction of 
the E- vector is orthogonal to the magnetic field). 
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The magn etic fields inferred from these observa- 
tions (e.g., [W erner et alJ[l988|; iMorris et aLlll992l : 
Hildebrand et al. 1993) run parahel to the Galac- 
tic plane. Recent linear polarization observations 
of the 450-/im continuum suggest a large-scale 
toroidal magneticjieid extending over a region of 
17 X 30 pc iNovak et al.ll2QQ3[ ). 

IChuss et all (|200 3^ found an interesting de- 
pendence of the magnetic field direction on the 
sub- mm flux. In low density regions, the field 
aligns generally perpendicular to the Galactic 
plane, while in high density regions, the field has a 
toroidal geometry. One explanation for this is that 
the global magnetic field in the early Galaxy was 
initially in a poloidal configuration, however the 
gravitational energy density in dense molecular 
clouds was strong enough to distort the poloidal 
field into a toroidal one, though it was insufficient 
in the low density regions. 

A model which can connect the poloidal 
and toroidal magnetic fields was proposed by 
Uchida et al.l ( 1985), and was exten d ed to a more 
realistic case bv 'Sh ibata fc Uchidal ( 1986h . This 
magnetohydrodynamic model w as developed to 
explai n the GC lobes found by Handal 
(|l984l ). Since magnetic flux is generally frozen 
into matter, differential rotation and infall can 
shear an initially poloidal field into a toroidal one. 
Consequently, a toroidal field is developed close 
to the Galactic plane, while the field is vertical at 
high Galactic latitude. 



The lUchida et al 



(|l985l ) model predicts that 
the toroidal component will generally be more 
dominant closer to the Galactic plane. A map of 
the direction of the magnetic field could therefore 
be a simple indicator of how well this model works 
in the GC. However, dust emissivity is high in 
dense, warm clouds, and thus observations of dust 
emission in FIR/sub-mm wavelengths are strongly 
limited to such regions, which show patchy distri- 
butions in the GC region. 

In this paper, we present near-infrared (NIR) 
polarimetry of point sources toward the GC cov- 
ering a much larger region of the sky than previous 
observations of this region. We demonstrate that 
NIR polarization can provide information on the 
magnetic field structure not only in the Galactic 
disk, but also in the central region of our Galaxy. 



2. OBSERVATIONS AND DATA RE- 
DUCTION 

We conducted NIR polarimetric observations 
of the GC with the SIRPOL camera on the night 
of 4 July 2006. SIRPOL consists of a single-beam 
polarimeter (a half-wave pl ate rotator un i t and 
a fixed wire-grid polarizer; iKandori et al.l l2006l ) 
and the NIR imaging camera SIRIUS (Simul- 
taneous Infrared Imag e r for Unbiased Survey 



[Nagashima et al. 19991 : Nagavama et al 



burvey ; 

Eooi, 



and is attached to the 1.4-m telescope IRSF (In- 
frared Survey Facility). SIRPOL provides images 
of a 7'.7 X 7'.7 area of sky in three NIR wavebands, 
J (1.25/im), H (1.63/im), and Ks (2.14/im), si- 
multaneously. The detectors are three 1024 x 
1024 HgCdTe arrays, with a scale of 0V45 pixel" ^ 
The filter system of IRSF /SIRPOL is similar to 
the MKO system (|Tokunaga et al.ll2002h . 

We observed a 20' x 20' area (nine SIR- 
POL fields) centered at the position of Sgr A* 
(17h45m4Q Qs^ -29°00'28:'0; J2000.0). We took 
10-s exposures each at 4 wave plate angles (0?0, 
22?5, 45?0, 67?5) at each of 10 dithered positions. 
The weather condition was photometric, with a 
seeing of ^V!2 (J), -I'.'l [H), and -I'.'O [Ks). 
The IRAF (Image Reduction and Analysis Facil- 
ity jil software package was used to perform dark- 
and flat-field corrections, followed by sky back- 
ground estimation and subtraction. 

3. DATA ANALYSIS AND RESULTS 
3.1. Polarization of Point Sources 

To obtain the photometric magnitudes and er- 
rors in the three bands, we used t he DAOFIND 
task in the DAOPHOT package (|Stetsonl Il987h 
to identify point sources in Stokes / images [/ = 
(/qO -h /22O5 + /450 + Iq^o^)/2]. The sources were 
then input to the ALLSTAR task for PSF-fitting 
photometry. About 10 sources were used to con- 
struct the PSF in each image. Each Stokes / im- 
age was calibrated with the photometric image of 
the same pos ition obtained in previo us imaging 



observations ( Nishivama et al. _ 2006ah. in which 



the standard star #9172 ( Persson et al.|[l998l ) was 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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used for calibration. We assumed that #9172 
has magnitudes of J = 12.48, H = 12.12, and 
Ks = 12.03 in the IRSF/SIRIUS system. 

In Fig. [U we show the H — Ks histogram (top 
panel) and the JHKs color-color diagram (bot- 
tom panel) for stars with photometric errors of 
less than 0.11 mag. Also plotted are loci of unred- 
dened giants and dwarfs. The arrow is parallel to 
the reddening vec tor, and its length corre sponds 
to Aks = 1 vciSig (jNishivama et al ] l2006al ). Con- 
sidering large extinction toward the GC, stars with 
small H — Ks should be attributed to foreground 
(disk) stars, and stars consisting of the strongest 
peak in the H — Ks histogram is attributed to ones 
in the Galactic bulge. 
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Fig. 1. — H — Ks histogram {top) and JHKs 
color-color diagram {bottom) for point sources 
with J, and Ks photometric errors of less than 
0.11 mag. The thick and thin curves re present the 
loci o f giants and dwarfs, respectively (iTokunagal 
2000l ). The arr ow indicates the Ak^ = 1 mag red- 



dening vector ( Nishiyama et al.]l2006ah . 



Astrometric calibration was performed, field 
by field, with reference to the positions of point_ 



used for the calibration. As a result of this as- 
trometric calibration, we obtained an rms of the 
positional difference of better than OVl for sources 
with a < 0.11 mag photometric error. 

The Stokes parameters /, Q, and U for point 
sources were determined from aperture polarime- 
try of combined images as follows. DAOFIND and 
APPHOT tasks were used for the point-source 
identification and the aperture photometry. We 
then obtained the intensity for each wave plate 
angle (/qQ, ^22° 5' ^45°' ^er^s)* ^^^^^ aperture pho- 
tometry with a small aperture radius gives a better 
photometric result than PSF-fitting photometry, 
aperture photometry was applied in the following 
procedure. The size of the PSF is slightly differ- 
ent among the four images due to variations of 
the seeing, hence we used different apertures for 
each image. The aperture diameters were set to 
be equal to 2 x FWHMs of the best fit Gaussian 
profile (GFWHM) determined by the PSFMEA- 
SURE task. The means of the adopted aperture 
sizes were -^.'2 (J), -I'.'l {H), and -I'.'O {Ks). 
The position angle off set of SIRPOL a was esti- 
mated to be a = 105° (|Kandori et all boOG) which 
was set to the origin of the angles 0°, 22? 5, 45°, 
and 67? 5. Based on the intensities of the four an- 
gles, we calculated the total intensity / and two 
"raw" Stokes parameters Q' and U' as 



(V +^22?5 



"^45° +^67?5)/2. 



^22:5 ^67:5- 
The polarization degree P and the position angle 
were derived by 



P = ^{Q'^^U'^)/I, 



= - arctan([/ / Q) , 



where 



Q = cos(2a) - sin(2a), U = sm{2a) + cos(2a). (1) 

The debiased P was finally derived by Pdb = 
\J P^ — bP^ where bP is the error of P, given by 



^P 



sourc es in the 2MASS point source catalog (Skrut skie et al 
I2OO6I ). Sources with photometric errors of less 
than 0.05 mag in 2MASS and our catalog were 




(2) 



We regard sources for which (P^ — (5P^) < as 
non-polarized source, and we do not consider such 
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sources further in this paper. The typical magni- 
tudes for 6P = 1% are 14.5 (J), 13.5 {H), and 
12.0 (Ks). 

Fig. [2] plots the degree and position angle of 
stars with a degree of polarization determined to 
have an accuracy better than 1%. The orientation 
of each bar gives the inferred direction of polariza- 
tion and the length of the bar is proportional to 
the degree of polarization. The coordinate offsets 
{') were measured with respect to the location of 
Sgr A*. 

Histograms of debiased P and in each band 
are shown in Fig. [3l The mean degree and angle 
are 6.3 % and 8?0 in the J band, 6.2 % and 13?6 
in the H band, and 4.3 % an d 14?4 in the Ks 
band. iKobavashi et al. detected i^-band 

emission of unresolved point sources at the central 
7' X 7' region and obtained an average degree of 
polarization of 5%, with position angles of 10° to 
15°, showing a good agreement with our results in 
the Ks band. 

The correlation between in the Ks band and 
H — Ks color is shown in Fig. [H To clarify 
the dependence of ^ on — Ks, we divided the 
H — Ks data set in bins of equal size (0.5 mag) 
and calculated the mean and the standard devia- 
tion of in each bin, represented by red crosses 
in Fig. m We can see a change of the mean po- 
sition angle at — Ks ~ 1.0: the mean angle is 
- 5° at < - < 1.0, while it is - 15° at 
H — Ks > 1 .0. Such a change has already been 
indicated bv IKobavashi et "aL ( 1983 ). and is also 
in good agreement with their results. 

We have identified two distinct populations in 
Fig. [3l stars with small P and small (typically 
^ ^ ^ stars with larger P 

and > 10° (see also Fig. |4]). It is conceivable 
that the former are nearby stars, and the latter are 
stars distributed in the Galactic bulge. The stars 
that corresponds to the strong peak in the Pj his- 
togram at ~ 2% have a color of H — Ks ~ 0.1—0.2 
(see also Fig. [T]) and Oj ^ 0° and thus are most 
likely to be nearby dwarfs. This peak can be found 
in the H band around Ph ~ 1%, but becomes in- 
visible in the Pks histogram. In the Ks band, 
most of the stars detected are red giants in the 
Galactic bulge, which constitute a strong peak in 
the Pks histogram. The strong peak at Pks ^ 4% 
corresponds to those at ~ 7% in Ph, and at ~ 10% 
in Pj, because this change of degree of polarization 



can be explained by the po wer law Pa cx: o f 
the interstellar polarization ('Nagata et al.lll993 ). 
The distinct populations, and the wavelength de- 
pendence of the polarization will be discussed in 
another paper (H. Hatano et al., in preparation). 

Since the vector maps in Fig. [2] are crowded 
and thus almost illegible, we show the i^s'-band 
mean vector map in Fig. [5l The mean degree 
and position angles were calculated using stars in 
a O.'S X 0.'8 grid with Pks/SPrs > 3.0. The vectors 
are superposed on the three color (J, Ks) com- 
posite image of the same region. At a first glance, 
most of the vectors are in order, which can also be 
seen in the Oks histogram, and are nearly paral- 
lel to the Galactic plane. Moving north-eastward 
across the image, the position angles slightly ro- 
tate clockwise. At a few positions where the num- 
ber density of stars is small, and hence strong fore- 
ground extinction exists, the vectors have irreg- 
ular directions, particularly at the northwestern 
corner. These irregularities might be explained by 
the inherent magnetic field configuration in fore- 
ground dark clouds. 

3.2. Separating Foreground Polarization 
and Galactic Center Component 

Polarization vectors of stars trace the plane-of- 
the-sky projection of the magnetic field, and po- 
larimetric measurements of stars of different dis- 
tances reveal the three-dimensional distribution of 
the magnetic field orientations. In our observed 
fields, we can detect stars in the Galactic disk 
and bulge. From the stars at the close side in 
the Galactic bulge (referred to hereafter as "blue 
stars" due to their relatively small reddening), we 
can obtain the degree of polarization and position 
angle, which are affected mainly by interstellar 
dust in the Galactic disk. The light from stars 
at the far side in the bulge (hereafter "red stars" ) 
is transmitted through the dust in the disk and 
the bulge. Therefore, using both blue and red 
stars, we can obtain the bulge (GC) component 
of the polarization. The procedure is as follows. 
(In the following process, only i^T^-band polariza- 
tion of stars is used.) 

As a first step, we divided the field into 10 x 10 
sub-fields of 2' x 2' and drew H — Ks histograms 
for each sub-field with stars of SH < 0.11, dKs < 
0.11, and H < 15.0, which are much brighter 
than the lOa limiting magnitude in the H band 
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Fig. 2. — Polarization of the Galactic center for 
stars with P > 0% and 5P < 1%. 2243, 7963, 
and 9661 stars in the J (top), H (middle), and Ks 
(bottom) bands, respectively, are plotted. The co- 
ordinate offsets (') were measured with respect to 
the location of Sgr A*. Each bar is drawn paral- 
lel to the E-vector of the measured polarization. 
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Fig. 3. — Histograms of degree of polarization 
(left) and position angle (right) in the J (top), H 
(middle), and Ks (bottom) bands, for stars with 
P > 0% and SP < 3%. 5795 (J), 17356 (iJ), 
and 18632 (Ks) stars are employed in these his- 
tograms. 



(16.7 mag). One of the H — Ks histograms is 
shown in Fig. [H upper left panel. Using the his- 
tograms, we evaluated a peak value of the his- 
togram {H — i^s')peak each sub-field. Fig. [71 
is an H Ys. H — Ks color-magnitude diagram 
(CMD) for stars with ^i^ < 0.11 and SKs < 0.11. 
This CMD shows that the criterion H < 15.0 is 
bright enough to avoid the influence of the lim- 
iting magnitudes {H ^ 16.7, Ks ^ 15.5) on the 
determination of {H — i^s')peak- 

Using the H — Ks color, we divided the stars 
with 6P < 3% into three sub-groups: "nearby" 
stars and "blue" and "red" stars in the bulge. 
We assume that nearby stars have a color of 
H - Ks < 1.0, because at - i^^ ^ 1.0 the 
number of stars drops and approaches a mini- 
mum (Fig. [H top panel) and a clear change 
of position angles can be determined (Fig. [4]). 
The "blue" stars are redder than H — Ks = 1.0 



H-Ks [mag] 



Riglit ascension (JSOQOjO) 



Fig. 4. — Position angle in the Ks band vs. 
H — Ks diagram for stars with SH < 0.11, SKs < 
0.11, > 0% and 6Pks < 3%. The red crosses 
represent the mean and standard deviations of 
in 0.5 mag width bins. 



and bluer than {H — ^5) peak- "^^^ stars with 
H — Ks > {H — ^5) peak selected as "red" 
stars (lower left panel in Fig. [6]). The average 
and standard deviation of {H — i^^jpeak 
sub-fields are 1.85 and 0.24, respectively; Fig. [8] 
shows the H — Ks histogram (top panel) of the 
red and blue stars, and their location in the Ks 
vs. H-Ks CMD (bottom panel). The blue stars 
have a peak at — Ks 1.5, while the red stars 
have a peak at — Ks ~ 2.0. The red and blue 
stars of all sub-fields are plotted, so their distribu- 
tion is overlapped in the CMD of Fig. [HI There 
are ~ 60 to ~ 340 stars with 5P better than 3% 
in each sub-field in the Ks band. 

As a second step to subtract the polarization 
originating in the bulge, Q' /I and U' /I histograms 
in the Ks band were constructed for the blue and 
red stars in each sub-field (upper and lower right 
panels in Fig. [6]). We calculated their means 
as < Q' /I >B, < Q' /I >R, < U' /I >B, and 
< U' /I >R. We then obtained the degree of po- 



Fig. 5. — iC^-band polarization vector map su- 
perposed on the three color (J, Ks) composite 
image of the Galactic center. The Galactic center 
is the bright yellow blob in the center. The mean 
Pks and Oks are calculated for each O.'S x 0.'8 
grid. Note here that stars with Pkq > 0% and 
Pks/^Pks > 3.0 are used for the calculation, so 
that the mean Pks might be overestimated. 



larization and position angle for the blue stars. 




and for the red stars. 




where < U/I >b and < Q/I >b, and < U/I >r 
and < Q/I >r are the Stokes parameters whose 
position-angle offset are corrected [see equations 
([1])]. We obtained P and for the red minus 
blue componen ts using the following equations 
([Goodrich! 1 19861 ): 



Or-b - arctan 




where (< U/I >r - < U/I >b) and (< Q/I >r 
— < Q/I >b) are also the Stokes parameters 
whose position-angle offset are corrected [see equa- 
tions dl])]. The errors of < Q'/I > and < U^I > 
were calculated from the standard error on the 
mean a / \/N of the Q' /I and U' /I histograms, 
where cr is the standard deviation and N is the 
number of stars. The error of Pr-b was then cal- 
culated by propagation of errors [see equation ([2])]. 



60 

50 

40 

30 

w 20 
(5 

o 
E 



1 


,|,,,,|,,,,|,,,,|,,, ,|, 

1 ; 

_n 


; } 






H-Ks [mag] 

Fig. 6. — Upper left: H — Ks histogram of a sub- 
field {l,h = with stars with SH < 0.11, 
SKs < 0.11, and H < 15.0. The arrow rep- 
resents the peak of the histogram. Lower left: 
H — Ks histograms of the same sub-field for blue 
(hatched) and red (dotted) stars with Pks > 0% 
and 6Pks < 3%. Upper and lower right: Q' /I and 
U'/I histograms for the blue (hatched) and red 
(dotted) stars shown in the lower left histogram. 



We show vector maps for Pb and Ob (blue bars) 
and Pr and Or (red bars) in Fig. [9l and the same 
map for Pr-b and ^r-b in Fig. [TOl The average 
of Pb and Ob are 3.8 % and 15? 1, and of Pr and Or 
are 4.3 % and 15?0, respectively (Fig. [TTl hatched 
and dotted histograms). Those of Pr-b and Or-b 
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Fig. 7. — vs. — iiTs' color-magnitude diagram 
for stars with SH < 0.11 and 5Ks < 0.11. The 
dashed line represents H = 15.0. 



are also obtained as 0.85 % and 16?0 (white his- 
togram in Fig. [11]), only for grids where the po- 
larization is detected with Pr_b/^Pr-b > 2. The 
averages for Ob and Or-b are similar, but their dis- 
persions are different: the standard deviation of Ob 
is 6?0, while that of ^r-b is 21? 5 (the average of 
^^R-B is 7? 6.) The small dispersion of ^b suggests 
that the long axis of the interstellar dust grains in 
the Galactic disk is well aligned perpendicular to 
the Galactic plane. The histogram of ^r-b has a 
peak at ~ 20°, which also roughly coincides with 
the angle of the Galactic plane. A similar result is 
also obtained for the i^-band polarization. 

In Fig. [TOl our Pr-b and Or-b are plotted 
overlaid on the P -vectors derived from FIR/sub- 
mm observations ( 



2QQ0l : IChuss et al 



Dotson et aT]l2000l : [Novak et al 
2003[ ). Although the observa- 



tions are restricted to the position of dense molec- 
ular clouds, we find good agreement between the 
FIR/sub-mm and NIK observations in spite of the 
difference of wavelengths and methods of deriving 
polarization. 
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Fig. 8. — H — Ks histograms (top) and Ks vs. 
H — Ks color-magnitude diagram (bottom) for red 
and blue stars with Pks > 0% ^ind 5Pks < 3%. 
Theoretical isochrones for different extinctions of 
^Ks = 2, 3, 4, and 5 mag, are shown by solid 
curves. The arro w indicates the Akc, = 1 mag 
reddening vector ( Nishivama et al.1l2QQ6a| ). 



4. DISCUSSION 

4.1. Previous Infrared Polarimetry toward 
the GC 

NIR polarimetry for diffuse emission and point 
sources tow ard the GC has been conducted since 
the 1970's (|Maihara et al.lll973h . From the po- 
larization angle aligned nearly along the Galactic 
plane, and the wavelength depend ence of polar- 



izatio n well fitted by a power law (|Nagata et al 
it has been interpreted that the polariza- 
tion is of interstellar origin do minated by dust that 



lies in the Galactic disk (e.g., Hough et al.l 1 19781 : 



Kobavashi et al.lll98Ql : ILebofskv et al.lll982h . 
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Fig. 9. — i^s'-band polarization map derived from 
the blue-star component (Pb & • ^lue bars) and 
red-star component (Pr & 6>r : red bars). 



One of the deepest NIR polarimetric ob- 
servations toward the GC was carried out by 
lEckart et al.1 (|l995h for a smah field of - 13'' x 13''. 
They established significant polarization for 160 
sources fainter than 13 mag in the K band. The 
mean flux- weighted polarization is 4 % at 25°, 
nearly parallel to the Galactic plane. A similar re- 
sult, 4.1 zb 0.6% at 30 ± 10°, wa s also o btained by 
lOtt et al.l (|l999h . lEckart et all (|l995h concluded 
that most of the polarizations are caused via ab- 
sorption by aligned dust grains in the Galactic 
plane. 

A change of the magnetic field configuration 
along the line of sight toward the GC has been 
pointed out by Kobayashi et aL (1983). The dia- 
gram of position angles in the K band versus H 



K fo r 15 discrete sources (Fig. 5. lKobavashi et al 



[1983) showed that the position angles are smaller 
and less ordered for sources with H — K < 1.0, 
while those with H — K > 1.0 are confined to 
a relatively narrow range around 20°. They con- 
cluded that this may indicate a change of the mag- 
netic field direction at a distance corresponding to 
H — K = 1.0 (5 kpc or more according to their 
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Fig. 10. — Ks-hdiiid polarization map derived 
from the Galactic center component (Pr-b & 
fe-B, Ted bars), where the polarization is detected 
with Pr_b/^Pr-b ^ 2. The polarization map at 
the center of our Galaxy derived from FIR/sub- 
mm observations {blue bars) is also shown. The 
length of the bars is proportional to the mea- 
sured degree of polarization, and their orienta- 
tion is drawn parallel to the inferred magnetic 
field direction. The FIR/sub-mm wavelengths 
da ta sets are froni 60 a m & 100 jam polarimetry 
by [Dot son et al. ' Nod), a nd 350 /xm p olarimetry 
by [Novak et al. (200^ and lChuss eTall ([2003, see 
also their Fig. 1). Some prominent radio filaments 
are shown as heavy dark lines. 
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Fig. 11. — Top: Histograms of degrees of polar- 
ization for Pb (hatched), Pr (dotted), and Pr-b 
(white). Bottom: Histograms of position angles 
for Ob (hatched), ^r (dotted), and ^r-b (white). 



calculation). 

We can also identify a change of the position 
angle in the Ks bands at — Ks ~ 1.0 (see Fig. 
H]). As shown in the color-color diagram in Fig. [TJ 
most of the stars with H — Ks ^1-0 have the color 
of giants (i.e., their unreddened positions are on 
the locus of giants), and the strong peak in the H— 
Ks histogram (the top panel in Fig. [1]) suggests 
that they are distributed in the Galactic bulge. 
Hence, this change of the position angle indicates a 
transition of the magnetic field configuration along 
the line of sight, between the Galactic disk and 
bulge. 



4.2. "Red" and "Blue" Stars and Their 
Separation in the Line of Sight 

To discriminate between foreground disk stars 
and those in the Galactic bulge, the criterion H — 
Ks = 1.0 is applied in our analysis because most 
stars with H — Ks > 1.0 are attributed as giants 
in the Galactic bulge. The histogram of H — Ks 
and the color-color diagram of point sources with 
photometric errors of less than 0.1 mag in the three 
bands are shown in Fig. [TJ Nearby stars can be 
found around the giants' and dwarfs' loci (thick 
and thin curves) without reddening, while most of 
the redder stars, which are located at — Ks ^ 
1.0 and J — H > 2.5, have a color expected for 



reddened giants. Toward the GC, the number of 
stars in the bulge is larger than in the exponentia l 



disk by a factor of ~ 50 ( Wainscoat et al J 1 19921 ). 



and thus most of the stars at — Kg ^1-0 and 
J — H > 2.5 are giants in the Galactic bulge. 

The bottom panel of Fig. [8] shows the lo- 
cation of the red and blue stars in the Ks vs. 
H — Ks color-magnitude diagra m. Th e theoret- 
ical Padova isochrones ( Gir ardi et al.l [2002) for 
solar metallicity with an age of 10 Gyr are also 
plotted in the color-magnitude diagram. The 
isochrones are put at the distanc e of the GC 
(7.5 kpc; iNishivama et al. 1 l2006bh and shifted 
along the re ddening vector {Akq / Eh-Ks = 1-44; 
INishivama et al. 2006a) with Ks extinctions of 0, 
1, 2, 3, 4, and 5 mag. The distances to the red 
and blue stars are difficult to estimate accurately, 
since we do not know the distribution of the inter- 
stellar dust along the line of sight; however, most 
of the giants in the Galactic bulge can be detected 
for at least Ks < 15, and for these stars, we 
can find a clear peak in the H — Ks histograms. 
Fig. [71 clearly shows that our observations are 
essentially equally sensitive to stars on the near 
{H -Ks < 2.0) and far {H - Ks > 2.0) sides of 
the Galactic bulge for H < 15.0. The clear peaks 
found in the H — Ks histograms are thus due 
to the spatial distribution of stars, not due to a 
distance effect. This indicates that {H — Ks)^^^^. 
roughly corresponds to a real peak of the spatial 
distribution of giants along the line of sight. The 
color-magnitude diagram of Fig. [8] also tells us 
that a large number of blue stars are distributed 
around the distance where 1.0 < Aks ^2.0, and 
most of the red stars are further than the distance 
which suffers from the extinction of Aks = 2.0. 

From the distribution of H — Ks color, we 
estimate the depth of the region where we have 
mapped the magnetic configuration shown in Fig. 
flOl For each sub-field, we calculated H — Ks 
color differences between the peak in the H — Ks 
histogram and the mean colors of blue and red 
stars, i.e., {H - Ks)peak - {{H - Ks)hiue) and 
{{H - Ks)red) -{H- i^5)peak. Thcsc color dif- 
ferences and corresponding extinction ^Ks^iue ~ 
1.44 X [{H - i^5)peak - {{H - Ks)une)] "and 

Aks^^, = 1.44 X [{{H - Ks)red) -{H- i^5)peak], 

where 1.44 comes from Aks / ^h-Ks = 1-44 
( Nishivama et al.l l2006at ). show the amount of 
dust extinction in the region where we have 



obtained polarimetric information. First, the 
amount of extinction toward the peaks in the 
H — Ks histograms can be calculated to be 
Aks,^^^ = 1.44 X [{H - Ks)peak -{H- Ks)o], 
where the mean intrinsic color {H — Ks)o of red 
giants is assumed to be ~ 0.2 (see Fig. [8]). This 
corresponds to the amount of extinction up to 
the GC in each sub-field. Next, the extinctions 
A^sbiue ^^sred couvcrtcd to thc actual 
distances from the GC, using a simple model. To 
m ake conservat i ve est imates, we use the model 
bv IPavies et al.l ( 1997 ). who derived a rather ex - 
tended distribution of dust. IPavies et al.l ( 1997 ) 
showed that the extinction to the GC associated 
with cool diffuse interstellar dust can be calculated 



usmg 



Cx 



Rq 



(3) 



where ~ 5.3 kpc is a scale- length of the dust 
distribution in the radial direction, (7 is a con- 
stant depending on the wavelength and dust den- 
sity, and Ro is the dista nce between the GC an d 
the Sun {Rq = 7.5 kpc: INishivama et al.ll2006bl ). 
The distances from the GC Xbiue and x^ed corre- 
sponding to ^Ksbiue ^^Sred' Tcspectively, can 
be estimated with the equations 



r/ad^^' 




We obtained the average and standard deviation 
of ^biue as 0.5 and 0.1 kpc, and those of Xred as 1.0 
and 0.5 kpc, respectively. We found a long-side tail 
in the x^ed histogram, which enlarges the average 
and the standard deviation. This estimation sug- 
gests that the polarization shown in Fig. [TOloccurs 
between the average distances of {Rq — 0.5) kpc 
and {Ro + 1.0) kpc from the Sun, arising probably 
from the central 1—2 kpc region of our Galaxy. 

In reality, the central part of the Galaxy har- 
bors a strong concentration of gas and dust called 
the "Central Molecular Zone" {R < 200 pc; 
[m orris fc Serabvnl Il996l ). which approximately 
corresponds to the concentration of star s called 
the "Nuclear Bulge" (R 5, 300 pc; Mez ger et al.l 
19961: ISerabvn MorrisI Il996h . According to 



Laun hardt et al. " (2002'), stars belonging to the 



Nuclear Bulge dominate in the central part of the 
Galaxy (see their Fig. 14). Therefore, a large 
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portion of the stars we have detected is located in 
the Nuclear Bulge, and the polarization at the GC 
originates mostly within the central few hundred 
pc. 

4.3. Magnetic Field Configuration at the 
GC 

As shown in Fig. [TOl we present the mag- 
netic field configuration at the central region of 
our Galaxy by discriminating between the polar- 
izations due to the disk and GC origin. The peak 
of the histogram of the position angle is ~ 20° 
(see Fig. [TT]), almost parallel to the Galactic 
plane. This coincidence, and the deficiency of 
^R-B around —60°, the angle perpendicular to the 
Galactic plane, show the basically toroidal geom- 
etry of the magnetic field. We cannot find a clear 
systematic dependence of position angle on Galac- 
tic latitude (from b ^ -0?27 to +0?18), suggesting 
that there is no systematic transition of the mag- 
netic field direction in this region. 

The direction of dust grain alignment at the 
GC has been investigated from polarized dust 
emission in the mid- and far-infrared and sub- 
mm wavelengths. The magnetic field implied by 
the emission is generally parallel t o the Galactic 
plane in the circuninuclear disk (e.g., Wern er et al.l 
19881 : iMorris et al.lll992l : iHildebrand et a l. 1993). 
As seen in Fig. [TOl the magnetic field configura- 
tion we obtained at the GC shows a go od agree- 



ment glo bally with tho s e obta ined by Do tson et al 




Novak et all (|2000f ). and iChuss et al 



which are the highest angular resolu- 
tion polarimetry data sets in the FIR/sub-mm 
wavelengths. The local features also show an 
excellent agreement: an X-shaped feature ex- 
tends from (Aa,AJ) = (— 3^ +100 down through 
(+5^00 (described in ^3.1.4. IChuss et all l2003h 
is also confirmed in our map. The configura- 
tion at M-0. 13-0.08 around (-1', -5') is also 
reproduced. The polarized FIR/sub-mm emission 
comes from molecular clouds, which are known 
to be located in the GC. Therefore we conclude 
that the position angles derived from our NIR po- 
larimetry represent the direction of the aligned 
dust grains in the GC. 

Chuss et al. I (|2003h has suggested that the mag- 



netic field aligns generally perpendicular to the 
Galactic plane in low density regions, while the 
field has a toroidal configuration in high density 



regions. They explain this correlation using the 
idea that in under dense regions, the magnetic field 
energy density can support itself against gravi- 
tational forces, preserving a primordial poloidal 
magnetic field. In over dense regions like molecu- 
lar clouds, on the other hand, the gravitational 
forces are strong enough to shear the magnetic 
field into a direction parallel to the Galactic plane. 
In this context, lower density regions should have 
a poloidal configuration. However, in our vec- 
tor map, the field shows a predominantly toroidal 
direction even at locations where the FIR/sub- 
mm emission is too weak for polarimetry [at the 
so utheastern c o rner i n Fig. [TOl see also Fig. 1 



Chuss et al.1 (|2003h and Fig. 2 in ^No vak et al 



(|2003l ).] The weak FIR/sub-mm emission suggests 
a paucity of dense clouds along the lines of sight, 
and the polarization in this direction can be con- 
sidered to be interstellar in origin. Hence the "in- 
terstellar" magnetic field at the GC is dominated 
by a toroidal configuration, and the primordial 
poloidal magnetic field, if it existed, is not pre- 
served today in this region. 

Several radio filaments exist in our observed 
field, and th ree of them are pr o minen t: the GC 



Radio Arc (|Yusef-Zadeh et al.l Il984t ). and the 



Northern and Southern Threa ds (also known as 
GO.08+0.15 and G359.96+0.09: lMorris Yusef-Zadehl 
1985), which are shown in Fig. [TOl Polarization 
studies have confirmed that the emission from the 
filaments is strongly linearly polarized, and that 
the internal magnetic field orientations are pa ral- 
lel to the long axes of t he filaments (Tsuboi et al.l 
1986[ iLang et al.lll999h . The simplest interpre- 
tation of these observations, combined with the 



disco v ery of more filaments 



20041 : lYusef-Zadeh et al.1 1200 



e.g., LaRos a et al.l 
is that poloidal 
magnetic fields are pervasive throughout the cen- 
tral few hundred pc. Although the polarizations 
originating from the filaments cannot be detected 
in our observation, those toward the surrounding 
fields can be detected. The filaments have a width 
10^' 



of less than ~ 10^' (|Lang et al.l 119991 ). and proba- 
bly have a depth similar to their width. On the 
other hand, the size of the grids of our analysis is 
2' X 2', and the polarization shown in Fig. [TOl is 
the average of ~1— 2 kpc from the GC along the 
line of sight. Hence most of the stars we detected 
do not show a polarization originating from the 
filaments. However, we can detect the average 



11 



polarization near the line of sight toward the fila- 
ments. As shown in Fig. [TOl most of the position 
angles of the grids including the filaments align 
nearly perpendicular to them rather than paral- 
lel. Therefore, Fig. [10] suggests that the average 
magnetic field has a toroidal configuration even 
around the sight-lines toward the filaments. We 
note again that the polarization is the average 
along the line of sight and does not originate from 
the area close to the filaments. 

4.4. NIR Polarization as a New Tool for 
Mapping the GC Magnetic Field 

We have shown that the polarization of starlight 
can be a prob e of the magnetic field near the GC. 
MorrisI (|l998l ) enumerated five different ways in 
which the magnetic field near the GC has been 
studied: morphology, polarization angle, Faraday 
rotation of the radio continuum, Zeeman effect, 
and polarized dust emission in FIR/sub-mm wave- 
lengths. The polarization of starlight has not been 
employed for such investigations. 

The optical polarization of starlight can trace 
the structure of the local magnetic field, but can- 
not detect stars near the GC due to large extinc- 
tion. NIR polarimetry has been carried out to- 
ward the GC prior to our observations (see §4.ip . 
but division of polarization into several compo- 
nents (e.g., originating from the Galactic disk and 
the central region) has not been done previously. 
The wide field-of-view of the NIR polarimeter SIR- 
POL, and the statistical treatment of tens of thou- 
sands of stars enable us to study the magnetic 
field near the GC; that is, the NIR polarimetry 
of starlight is a new way to investigate the mag- 
netic field. 

NIR polarimetry has the advantage of provid- 
ing information about the magnetic field at loca- 
tions where FIR/sub-mm emission is weak. NIR 
polarization of starlight is attributed to extinc- 
tion along the line of sight by aligned dust grains, 
while FIR/sub-mm polarization is due to emission 
from the aligned dust. Hence, NIR polarimetry 
can investigate the magnetic field in regions where 
FIR/sub-mm polarimetry is absent, if background 
stars exist. This advantage is clearly shown in Fig. 
[TOl We have detected polarization at positions 
where blue bars are not shown. The distribution 
of the position angles including such low emission 
regions shows a globally toroidal magnetic config- 



uration at the GC. 

5. SUMMARY 

We have measured the near-infrared polariza- 
tion of point sources toward the Galactic center 
(GC) in the 20' x 20' region centered at Sgr A*. 
The difference in the Stokes parameters between 
stars at the close and far sides of the GC re- 
veals the polarization originating from the central 
1—2 kpc region of our Galaxy. The distribution 
of the position angle for the central region shows 
good agreement with those obtained from polar- 
ized emission of dust in the GC, showing that the 
near-infrared polarization of point sources can be 
used as a tool to investigate the magnetic field 
configuration of the GC. The position angles have 
a peak at ~ 20°, which is almost parallel to the 
Galactic plane, suggesting a global toroidal mag- 
netic field in the region. 
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